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' Abstract. For a class of singular potentials, including the Coulomb potential (in three 

. . and less dimensions) and V(x) = g/x 2 with the coefficient g in a certain range (x being 

a space coordinate in one or more dimensions), the corresponding Schrodinger operator 
is not automatically self-adjoint on its natural domain. Such operators admit more than 
one sclf-adjoint domain, and the spectrum and all physical consequences depend seriously 
on the self-adjoint version chosen. The article discusses how the self-adjoint domains can 
be identified in terms of a boundary condition for the asymptotic behaviour of the wave 
functions around the singularity, and what physical differences emerge for different sclf- 

• adjoint versions of the Hamiltonian. The paper reviews and interprets known results, with 
the intention to provide a practical guide for all those interested in how to approach these 
ambiguous situations. 
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\Q ', 1 Introduction 

\o 

OO 

• Let us consider a quantum mechanical Schrodinger Hamiltonian 



where the real potential V is singular, for example, V(r) ~ - or V(r) ~ 4. To keep the 
discussion technically simple, let our configuration space be a one dimensional interval or half 
line or line. Remarkably, the principles and methods to come will be valid for higher dimen- 
sional configuration spaces, too, it is only the amount of technicalities that increases with the 
dimensions. It is also important to note that the half line case plays a direct role in a number of 
higher dimensional problems as well. For instance, in a higher dimensional setting with a central 
symmetric potential that is divergent at the centre, one can perform a separation of variables 
in a spherical coordinate system, and the ambiguity that will be our topic here appears in the 
radial part of the system, which is actually a half line problem. 

Now, in a straightforward manner, let us innocently search for all the eigenfunctions of 
this differential operator. Then, depending on our "luck", we may encounter some surprising 
difficulties Namely, we may find that, actually, there are "too many" eigenfunctions, in 
the sense that they are not mutually orthogonal, are not linearly independent, and form an 
overcomplete system rather than a basis in the Hilbert space of square integrable functions. 



''This paper is a contribution to the Proceedings of the 3-rd Microconference "Analytic and Algebraic Me- 



thods III". The full collection is available at http://www.emis.de/journals/SIGMA/Prague2007.html 



1 Naturally, this will not be a question of "luck", and we will see soon the mathematical criterion that tells 
that which potentials cause these difficulties. 
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Moreover, we can observe that there are so many eigenfunctions that even the eigenvalues are 
not restricted to real values but any nonreal A € C \ R also proves to be the eigenvalue of some 
(square integrable) eigensolution of the eigenvalue problem. Apparently, the self-adjointness of 
our, naively self-adjoint, Hamiltonian is seriously challenged. This can also be seen from that, 
for generic wave functions ip and x> applying integration by parts, 

(Hip, x) — (ip, Hx) = surface terms 7^ (2) 

will be observable. 

Now, if we don't have self-adjointness then we have no spectral theorem, no physical in- 
terpretation, and no unitary time evolution. Therefore, we are heavily motivated to restore 
self-adjointness, if only possible. 

This last observation ([2]) may suggest us that the problem is created at the boundaries, 
especially at the location of singularity. We can try to cure the situation by requiring that 
all wave functions should vanish or decrease fast enough as we approach the singularity so 
that the surface terms tend to zero. Unfortunately, we will find that this way we lose all the 
eigenfunctions. Hence, this requirement cannot help in solving the problem. It seems that we 
need to keep some of the eigenfunctions and to dispose the rest of them. 

Actually, the usual reaction to the problem - as many textbooks treat the three dimensional 
Coulomb problem, for example - is to declare that not all eigenfunctions are "acceptable" . Sud- 
denly, some extra ad hoc requirement is invented, which does not follow from the axioms/ (general 
principles) of quantum theory but is introduced on the fly. Unfortunately, it is beyond the scope 
of this discussion to analyze those various conditions in detail. To give only a summary of the 
result of such a careful and honest analysis, one can reveal that 

• those requirements are physically questionable, 

• they have mathematically limited availability, 

• different conditions can lead to different results, 

• many consistent and valuable quantum mechanical models are lost by those requirements. 

The situation is well illustrated, for example, by the paper p], which, for the attractive Coulomb 
potential in one dimension, gives a critical overview of various existing treatments and their 
various different results - ands adds another new approach, which is also questionable both 
mathematically and physically. 

However, we can follow another philosophy as well. Instead of the urge to choose, we can 
accept that the potential itself does not fix the model uniquely. When an application forces us 
to fix the ambiguity, some additional physical information will be needed. That information 
can come from experimental measurement or, if available, from some additional theoretical 
knowledge about the concrete case. 

Therefore, for these potentials, let us search for all the quantum mechanically allowed cases 
(self-adjoint Hamiltonians) related to our initial Hamiltonian. Then, when in a concrete problem 
we need to choose one, 

• knowing what possibilities exist may help to choose, 

• if we have some additional information then it will be easier to utilize it, 

• if measurement is needed to decide then we know what to measure and what to fit to the 
experimental data. 

To carry out finding all the cases, it is advisable to review first what mathematics knows 
and tells about self-adjointness and about the possible ambiguity in it. 
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2 Self-adjoint extensions of symmetric operators 

We will make use of the following definitions and theorems^ (the validity of which is, nevertheless, 
not restricted to differential operators). 

Let A be a linear operator that is defined on a dense subset T>(A) of a separable Hilbert 
space H. 

The adjoint A + of A is defined on those vectors x £ ^ for which there exists a^£W such 
that 

(Aip,x) = (Ax) for y^eV(A), 

and A + is defined on such a^as A + x ■= X- 
A is called symmetric if 

{AiP, X ) = (?P,Ax) for ViP, X £V{A). 

The adjoint of a symmetric A is always an extension of it (i.e., T>(A) is a subset of T>(A + ), and 
A and A + act the same way on D(A)). 

A is self-adjoint if A = A + . (Which includes that their domains coincide.) A is essentially 
self-adjoint if it admits a unique self-adjoint extension. 

Remarkably, not the symmetric but only the self-adjoint operators are those for which the 
spectral theorem holds, and which are in a one-to-one correspondence with the strongly contin- 
uous one-parameter unitary groups U(t) on 7i (to any U(t), there is a unique self-adjoint A such 
that U(t) = e~ tAt ). Both these properties play an important role in the physical interpretation 
of quantum mechanics. 

Next, from now on, let A be a symmetric operator that is closed (for symmetric operators, this 
simply means the requirement A ++ = A] any symmetric operator admits a closure - a minimal 
closed extension - which is actually nothing but its double adjoint). Let us also fix an arbitrary 
nonreal number A. 

The deficiency subspaces £\ and 8\* are defined as the eigensubspace of A + with respect to 
the eigenvalue A and the complex conjugate eigenvalue A*, respectively. The deficiency indices n\ 
and n\* are the dimension of £\, resp. £\*. A is self-adjoint if and only if n\ = n\* = 0. 

A admits self-adjoint extensions if and only if its deficiency indices are equall, n\ = n\* =: n. 
The self-adjoint extensions are in a one-to-one correspondence with the unitary maps from E\ 
to E\* . Each unitary C/n : £\ — > £\* characterizes a self-adjoint extension Ajj n as the restriction 
of A + to the domain 

V(A Un ) = {Vo + tpx + U^x I € T>(A), Va € £a}- (3) 

The unitary maps Un act on an n dimensional space so they can be parametrized by n 2 real 
parameter^. 

After the characterization ([3]) given by von Neumann, let us see a more recent alternative 
description of the possible self-adjoint extensions, the so-called boundary value space approach: 

If there exists an (auxiliary) Hilbert space 7i\y (necessarily n dimensional) and two linear 
maps r 1; T 2 : V(A+) -> 7i h such that, for VV>, X G V(A+) 1 

(A+ip, X ) - (</,, a + x) = (r^, r 2X ) b - (r 2 v, r lX ) b , (4) 



2 All mathematical ingredients quoted here and hereafter are taken from the sources [2] [3] [4] [5] [6] . 

3 For our differential operator Hamiltonian ([T]) with a real potential, this will always be the case, since such 
an operator is invariant under complex conjugation. 

4 For our differential operator Hamiltonian in one dimension, n is finite (unless the potential admits infinitely 
many singularities). In higher dimensions, n can be oo, which happens when the singular locations are not finitely 
many points but lie, say, along a line. 



4 



T. Fiilop 



and to any two 1 ^> 1 , ^ 2 £ TCb there exists a ip G 2?(^4 + ) satisfying 

r^ = *i, r 2 v> = * 2 , (5) 

then there is a one-to-one correspondence between the unitary maps U £ UiTi^) and the self- 
adjoint extensions of A, where a C/ describes the domain of the corresponding Ajj as 

v(Au) = {ipe v(A+) | (u - i b )r^ + i(c/ + i b )r 2 ^ = o}. 

Heuristically, the second condition, related to ([5]), says that the auxiliary boundary Hilbert 
space 7ib must be a smallest one among those fulfilling the first condition, @. If an 7i\> is 
suitable for ^ then some appropriate trivial "enlargement" , extension of it can also be suitable, 
e.g., by orthogonally adding some other Hilbert space to it, so this second condition is to ensure 
the efficiency of the description by removing the redundancy. 

One can find that the former characterization of self-adjoint domains can be considered as a 
special case of the latter, with Ti^ = E\. Another note to make is that, for a fixed TC^, the choice 
of appropriate Ts is not unique. Hence, it is not unique that which U provides which self-adjoint 
version. Therefore, one should not - at least at the general level - attribute any distinguished 
meaning to that, having one choice of TC^, T 1 and T 2 , which self-adjoint version is indexed by 
which unitary operator. (See more on it later.) A third remark is that, although the boundary 
value space approach may appear a rather abstract equipment at first sight, in applications we 
can find it simple, friendly and handy. Actually, historically it has been worked out to directly 
suit the special cases of differential operators where the "boundary values" ([5]) are indeed the 
limiting values of the wave function and its derivative at the boundary of the configuration space 
- or appropriate combinations of them. Soon we will see examples that show how practical this 
approach is for differential operators. 



3 Finding the self-adjoint Hamiltonians 




Applying the content of the previous section to our 
differential operator of the type — J^- A +V, we first 
need to specify an initial domain on which it is symmet- 
ric. We can choose, for example, those wave functions ip 
which admit a continuous second derivative, vanish in 
a neighbourhood of any singularity of the potential 
and near any finite endpoints of our configuration 
space (if applicable), have compact support (so that 

they also vanish "in a neighbourhood of infinity", for the case when the configuration space 
has +oo and/or -co as "endpoint"), and ip, Hip are square integrable. See the figure for an 
illustration. Such a domain is dense in L 2 . 

The closure of this symmetric operator will have a bit generalized domain, where the smooth- 
ness property is weaker (ip and its derivative must be absolutely continuous), and the wave func- 
tions do not need to vanish identically around any "problematic" location but only to decrease 
fast enough so that the limiting surface terms of ([2]) emerging at the boundaries and singularities 
will be zero for any pairs ip, x- At last, the adjoint of this closed operator will have the further 
generalized domain in which there is no restriction on the behaviour around the "problematic" 
locations. 

Now, to specify the self-adjoint domains, which lie in between the symmetric domain and the 
adjoint domain, let us use the boundary value space description because that requires the least 
calculational efforts, and because there are known recipes how to find suitable candidates for the 
needed ingredients H.^, T 1 , T 2 . In what follows, we will see some instructive examples how the 
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method works in practice. The first two examples are chosen to be simple to make the essence 
of the procedure apparent. In the meantime, they will already exhibit many of the general 
physical properties that typically arise when one has a singularity-induced or boundary-induced 
ambiguity in self-adjointness. 



4 Free particle on a half line 

As the first example, let us consider a free particle moving h i d2 

on a half line, which we wish to be bordered by a perfectly 
reflecting boundary to ensure the conservation of probability 
and, correspondingly, a self-adjoint domain for the Hamil- 
tonian. (See the figure for an illustration and for the notations.) 

The most well-known way to reach this is to impose the Dirichlet boundary condition ip(0) = 0. 
However, other boundary conditions also provide self-adjoint possibilities, and physically it is 
a natural idea that there can be various different types of reflecting walla^l Therefore, exploring 
all the self-adjoint versions means to find the variety of perfectly reflecting walls that can be 
accounted for in the scope of this quantum mechanical model. 

Starting with a symmetric domain and considering the corresponding adjoint domain as 
discussed in the previous section, integration by parts allows us to write, for the adjoint operator, 

where W[-,-] denotes the Wronskian. In principle, we could expect a surface term related to 
x — > oo as well (since it is not immediately clear whether the square integrability of tp and H + ijj 
plus the degree of smoothness ip has provides a fast enough decrease towards infinity to send 
that term to zero): it will be explained in Section [6l at a general level, why it is absent indeed. 

Comparing this formula with (pEJ) shows that we can comply with the boundary value space 
approach with the choices 

n h = c, = ^(o), r 2 v = l ^'(o). 

Here, we have introduced an auxiliary nonzero real length L , which does not play any principal 
role in the question of self-adjointness but is needed merely on dimensional grounds. Namely, the 
formalism requires T 1 and T 2 to have the same dimensions while ip and ip' differ in a dimension 
of length^. Also, we have dropped the factor — ^— (that is, (|4|) has actually been fulfilled for 

The other, "efficiency", condition needed for the boundary value space approach is also not 
hard to check. This point we are actually going to discuss later in Section [6l again at a general 
level. 

Having fulfilled the requirements of the boundary value space method, we are allowed to 
harvest the fruit: the possible self-adjoint domains form a one-parameter family, indexed by 
U £ li(V), in a parametrized form, U = e l ® , d G [0, 2tv), via the boundary condition 

(e™ - 1) V(0) + i(e" + l)L o i/(0) = 0. 

We can imagine, for example, a long and thin nanowire, with a small piece of some material attached to its 
ends. Then, near one such end, quantum effects in the range of wavelengths that are much larger then the width 
of the wire as well as the size of this attached blob but smaller then the length of the wire may be well described 
by one of the "free particle on the half line" self-adjoint Hamiltonians, presumably different ones for different 
materials attached. 

6 For a full theory of how physical dimensions can be formulated and treated mathematically - via one dimen- 
sional vector spaces and their tensorial products, quotients and powers - see [7] . 
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We can rewrite this condition in the more simplified form 



ip(0) + Lil/(p) = 0, 



■d 

L = L cot— G (—00, 00) U {00} 



as well. Hence, the quantum mechanically allowed reflecting walls are characterized by an arbi- 
trary real length parameter. The Dirichlet case corresponds to L = 0, the Neumann boundary 
condition ^'(0) = is the case L = 00, and the remaining cases, containing a nontrivial com- 
bination of the wave function and its derivative in the boundary condition, are often called the 
Robin boundary conditions. 

It is instructive to see how the spectral properties, and correspondingly all physical properties, 
depend on the boundary parameter L. Omitting technical details, solving the eigenvalue problem 
for a given boundary condition yields [9] 



E > : if k {x) 



E < : Abound (x) 



-e- x ' L , 



E 



,2i<5 fc 



1 - ikL 
1 + ikL 1 



bound 



2mL 2 



(only if L > 0). 



Therefore, for example the existence and energy of a boundary bound state is heavily in- 
dependent. The L-dependence of the scattering phase shift can also be visualized by, e.g., the 
time delay, which is the difference between the time 
when the peak of an incoming wave packet reaches the 
wall and the time when the peak of the reflected packet 
leaves the wall. For an incoming packet concentrated 
in wave number around the value —k, the result is [HJ[9] 



-2mL 



hk(l + k 2 L 2 ) 



reflected 
k 



1 \ 
1 \ 
1 \ 
/ \ 
/ \ 




which seriously depends on L (even its sign is decided by the sign of L). For further L-related 
effects and aspects, the Reader is asked to consult [8j E] - 



5 Free particle on a line with a point interaction 

In the next example, the particle can move on a line 

freely, except at one point where some object or ^ * x 

disturbance or short-range potential-like effect resides 
and can perform some nontrivial action on the particle. 

In this setting, the widely known example is the Dirac delta potential but the full family of 
possibilities again proves to be larger, and can be explored again via finding the possible self- 
adjoint versions of the initial Hamiltonian. 

We can proceed similarly to the half line case, writing 

(H+ip, x ) - ty, h + x ) = ~ {(VY - ^*'x)(+o) - (VV - r'x)(-o)} , 

2m 

and finding agreement with through 
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Correspondingly, the family of self-adjoint possibilities can be indexed by U £ U{2), via 

(u - lc^Trf + i(u + i c2 )r 2 v = 0. 

Similarly to L of the half line case, here, the four real parameters parametrizing the U s can be 
chosen as two length scales L + , L_ and two angles. Roughly speaking, the reason for two length 
parameters is that, now, the singular point has two sides. In parallel, one of the angles expresses 
the mixing between the left and the right side, and the other characterizes the strength of 
a short-range vector potential and thus causes a phase jump in the wave function at the singular 
poin10. In the case of zero mixing, which happens when the matrix U is diagonal, there is no 
probability flow crossing the singularity, the two half lines physically decouple, and the system 
becomes a sum of two independent subsystems, two half line systems indeed (one with L = L + 
and the other with L = L_). 

Just as the sign of L decides the number of bound states in the half line model, here the signs 
of L_|_ and L_ govern the number of bound states, which can thus be 0, 1 or 2. The scattering 
properties also depend on the boundary parameters, causing time delays, and making the point 
object to act as a low-pass or high-pass filter depending on the values of the parameters. See 
[3 [ini [HI CE2] for further physical effects caused by the properties of the pointlike object. 



6 Singular potentials 

Now, we enter the situation having a potential that possesses not an abrupt ( "hard" ) singularity 
like a Dirac delta but a continuously developing one like for the Coulomb potential ("soft" 
singularity). The complication with respect to the abrupt cases is that, here, the wave functions 
in the adjoint domain and their derivative are typically diverging when approaching the location 
of singularity. Therefore, we cannot choose such limits as boundary values. However, there exist 
some finite numbers hidden in the asymptotics of the wave functions, which can be extracted 
as follows. Let us first restrict ourselves to one side of the singularity. 

One can start with observing^] that, for any ip, \ hi the adjoint domain, W[ij}*, %] has a finite 
limit when approaching the singularity. Next, fixing an arbitrary real value Eq, let us consider 
two real eigenfunctions (p^ 1 ', of our Hamiltonian as a differential operator (i.e., omitting 
the requirements related to square integr ability) for the eigenvalue Eq, two such solutions that 
satisfy VF^ 1 ), (p^] = 1. We will call them reference modes. Two possibilities exist: 

1. Both reference modes are square integrable in a (one-sided) neighbourhood of the location 
of the singularity - this is called in the mathematical literature the limit- circle cas^ 
(example: £ (-■£, + ^) on (0,oo), E = 0, <pV>( X ) = x 5 / 4 , ^(x) = -fx" 1 /*), 

2. At most, only a specific linear combination of theni 1 "^! is square integrable around the 
singularity - this is called the limit-point case (example: ^ ( — ^ + ^x^~ ) > on (0j°°)> 
E = 0, <p<»(x) = x 7 / 4 , ^ 2 \x) = -fx" 3 / 4 ). 

It turns out that, if our problem is in the limit-point case then, for any two wave functions tp, % 
in the adjoint domain, tends to zero at the singularity. Hence, a limit-point singularity 

does not create nonzero surface terms in ([2]) and, consequently, does not induce an ambiguity 

7 When replacing the line with a circle, this fourth parameter corresponds to the magnetic flux driven through 
the circle. 

8 For the mathematical results quoted in this Section, see, e.g., [3j[4][5]. 

9 The name has nothing to do with the geometry of our configuration space but is because of some historical 
and technical reasons. 

10 Up to an overall complex factor, of course. 
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in self-adjointness. As a special case, this is the reason why, in the "free particle on a half line" 
system, we did not need to worry about the surface terms towards infinity. Note that an infinite 
"endpoint" is always to be treated as a singular point, which can be seen, for example, from that 
the change of variable x — ► - brings it into finite (to zero), and via the accompanying unitary 
transformation i/j(x) := -ip(-) our symmetric differential operator is mapped into a one that is 

singular at zero (e.g., — ^:pr is mapped to — (^jg) + 2a; 2 ] , and see the last paragraph 

of this Section). 

In the limit-circle case where the surface terms do not vanish, what can be shown is that any 
ip in the adjoint domain is asymptotically similar to an appropriate linear combination of the 
reference modes, in the sense that 

tP( x ) = [c^ + V ^(x)]ip^(x)+ [c^ +r ] { - 2 \x)]^ 2 \x), 

where rj^\x), r]^ 2 \x) vanish towards the singularity, and the limit numbers can be 

read off as the (always finite) limiting values of — VK^ 2 ), ip], respectively VF^ 1 ), ip], at the 
singularity. Furthermore, this asymptotic similarity proves to be strong enough for that, when 
we want to calculate the limit of a W[y>*,x] at the singularity, we can replace ip(x) in it with 
c ( 1 ) ( ^( 1 )( a; ) + c^(p^(x) (and x can also De replaced with its similar approximation). 
Adding 

w[ip*,x) = w[^\ii>]*wy?\x] - w[ip® ,i>]*w[<pQ-\x\, 

which is nothing but a simple identity about determinants (recall VFf^ 1 ) , (p^] = 1), we find 
that the surface terms can be expressed in the desired form in terms of finite quantities so these 
limit numbers can be used for the purposes of I\, T 2 - The examples coming in the following 
sections will show this in close detail. 

Now we can see the big practical advantage of the boundary value space method to using 
von Neumann's characterization directly. Namely, here we need to check square integrability 
of eigenfunctions only in a local neighbourhood of the singularity while, in the von Neumann 
approach, we need to do it over the whole configuration space. Besides, there are some additional 
benefits as well. 

First, in the von Neumann method, we need to find eigenfunctions for a nonreal eigenvalue, 
while here for a real one. For example, many eigenequations become simpler for the eigenvalue 
E = 0. 

Second, we actually need only one solution, as the condition f^] = 1 allows us to 

find another one in a form 

^*):=^W* — f^. 

Jx [ipW(x)\ 

Third, we need the reference modes only approximately, to only such a preciseness that the 
limit of their Wronskian with any ip can be determined, plus that their local square integrability 
can be checked. In practice, this usually means that we need to know their leading and first 
subleading asymptotic behaviour. 

In case our singularity is a so-called regular singular point then the asymptotic solutions are 
in fact known from the Frobenius method. It is usually advantageous to choose the regular 
solution for one of the reference modes. 

As a last virtue to mention, it is also apparent that the need for a boundary condition, as 
well as its form, is decided only locally around the singularity and not globally along the whole 
configuration space. Using the von Neumann approach with the deficiency eigenfunctions this 
point would also remain hidden. Actually, if we are on an interval bordered by two limit-circle 
singular endpoints then we are allowed to use different reference modes at the two endpoints. 



Singular Potentials in Quantum Mechanics 



9 



Concerning the efficiency condition not yet discussed, the mathematical literature ensures 
that this will also be fulfilled with the r i; T 2 chosen above. Indeed, it can be shown that, to 
any two complex numbers , c*- 2 -* , one can find a tp in the adjoint domain whose limit numbers 
are these and . Observing that a regular endpoint (where the potential does not diverge) 
behaves the same way as a limit-circle singular endpoint from all relevant aspects - we can use 
X^\x) = —x and x^ 2 H x ) = 1 as approximate reference modes - one can obtain that, in the 
special case of a regular endpoint, the limits of tp and ip' as boundary values do satisfy the 
efficiency condition. 

The presented description of the self-adjoint domains via reference modes contains some 
arbitrariness. One such freedom is in the value of Eq. Since the term containing Eq in the 
eigenequation is overwhelmed by the diverging potential term, it is plausible to expect and can 
actually be proven that the asymptotic behaviour of the reference modes that decides the limit 
numbers of the wave functions does not depend on Eq. The other arbitrariness is how the two 
independent reference modes are chosen for a given Eq. This latter uncertainty, which is an 
SL(2,M) amount of freedom, really influences the characterization. (As a special case of what 
has been mentioned about the non-uniqueness of Ts in Section [2j) The family of self-adjoint 
domains is nevertheless the same, what changes is only that which domain is indexed by which U. 
Correspondingly, parameters chosen to parametrize U E ^(2) may also change. Naturally, it 
is advantageous if we are able to introduce parameters in such a - reference mode-dependent - 
way that they remain the same under changing the reference modes. 

A last remark is that, should our Hamiltonian have a kinetic term with ^{p( x )^) instead 
of we only need to replace ^ with p(x)^ in the above formulas and considerations. Any 
possible zeros oi p(x) are also to be considered singularities which, with this replacement, can be 
treated analogously to the singularities of an operator with constant p. An example is provided 
by the operator — [-r- (x^-^) + 2x 2 ] mentioned above. The zero of p(x) = x 4 , x = 0, indicates 
a singular point there. The reference modes are to be normalized to pVFf^ 1 ), f^] = 1, and can 
be chosen for Eq = as ip^\x) = x~ l , ip ( - 2 \x) = —x~ 2 . These reference modes indeed behave 
singularly at x = 0, and show that, in this case, the singularity is of limit-point type as neither 
of them is square integrable in a finite neighbourhood of x = 0. 



7 Singular potential on a half line 



Let us again take two examples - actually, the generalizations of our previous two examples 
to apply the knowledge collected in the previous section. 

First, let us consider a half line again, but now with a 
potential that diverges in a limit-circle way at the finite 
endpoint. Based on the findings and notations introduced 
in Sections 2] and [61 we do not need much explanation why 
the boundary value space approach and the corresponding 
characterization of self-adjoint domains can be given aa 11 ! 



1 = 



n h = c, 



r 1 ,/j = w[^ 1 \^}(+o), r 2 v = w[^ 2 \ip](+oy, t^ + lt 2 ^ = q. (7) 



To be concrete, we may choose the three dimensional Coulomb problem, with the Hamiltonian 



H 



2m 



A+ : - 



r 



11 An Lo may be needed to introduce here or there, because of the already mentioned dimensional reason, but 
it is also possible that a coefficient from the potential can take this role. L in @ may also have a dimension 
different than length, depending on the conventions we choose. 
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After the usual separation of variables in spherical coordinates (including the radial mapping 
L 2 ((0, oo), r 2 dr) — > L 2 ((0, oo), dr)), the radial part, 



2m 



d 2 | + | g 

^^^■2 ^»2 



rad \ 



will have an ambiguity in self-adjointness in the I = angular momentum channel because the 
centre is a limit-point singularity for the radial Hamiltonians H^ l X for / > but is a limit-circle 
one for / = 0. 

For I = 0, a satisfactory approximation for reference modes near r = is 



(pP'(r) 1 + grln|<7|r. 



(8) 



These can be used in the boundary condition ([7|). For L = 0, the widely known and usually 
considered case yields, in which the ips allowed by the boundary condition are regular at the 
centre. In the other cases L / 0, a singular component is present in ip so ip' diverges at the 
centre. Experiments make us quite confident in that, when the potential is intended to express 
a purely electromagnetic interaction, we should choose L = in addition. However, for a physical 
situation where some additional, short-range, interaction is expected to be present between the 
centre and the particle, a model with I / is a good candidate. Mesic atoms, in which 
not electrons but some negatively charged mesons like tt~ move around the proton-containing 
centre, provide an example for this, because of the short-range strong force also acting between 
the protons and mesons. 

For generic L, the condition for bound states can be 
found [13] to be 



F(f) 



gF 



IsJ-lmEjU 2 



1 

I 



(9) 



with 



F(£) = + i) - In |e| - Ti - *(1) - * (2), 

where VI/ denotes the di-Gamma function. 
If our potential is attractive, g < 0, then, for L = 0, we obtain the well-known bound state 
energies 



F - R 

&n — n, 



n= 1,2,... 



(L = 0). 



However, taking for example the case L = oo, the solutions of the transcendental equation ([9|) 
give a different infinite sequence of bound states in the I = angular momentum channel: 



E,, 



R 



(n - c n f 



n = 1, 2, . . . (L = oo) 



with 



ci « 0.5130, c 2 w 0.4879, 



c 3 « 0.4857, 



0.4844 (limit). 



The scattering properties - like the scattering length and phase shift quantities - and other 
physical aspects are similarly considerably influenced by the value of the self-adjointness para- 
meter L. 
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8 Singular potential on a line 

Next, let us again extend our configuration space to a line, 
on which the potential admits a singularity that is of limit- 
circle type from both directions. Similarly to the "free line 
with a pointlike singular object" example seen in Section^ 
the two sides of the singularity double the dimension of the 

boundary value space, and we also need to double (to extend) the reference modes. Having 
seen © and (J7|), no wonder that we choose 

and obtain the boundary conditions, with U £ ^(2), as 

(u - i C 2)r^ + i(u + i c2 )r 2 v> = o. 

Here again, a certain subfamily of the family of self-adjoint Hamiltonians contains the separated 
cases where the system decouples to two independent half line subsystems. This means that, in 
those cases, the particle bounces back from the singularity. This phenomenon is present for at- 
tractive potentials as well, which is surprising for the physical intuition. It creates a feeling that, 
in those cases, something repulsive must happen at the location where the potential diverges. 

On the other side, in the nonseparated cases we can 
observe reflection-transition, time delay and filter phe- 
nomena similarly to the case of a pointlike interaction on 
a line. This includes that we find tunneling through a re- 
pulsive potential as well [H] , which is much more surprising 




for such a "fat" infinite potential barrier than for a "thin" pointlike object. Again, one may 
try to illustrate these cases physically by that some attractive effect is present at the location 
where the potential is undefined. 

Many other physical properties, like integrability for Calogero-type potentials whose singular 
term has g < |^ |15j . are also dependent on the self-adjoint version choserl^l. These depen- 
dences might be utilized in the future to build quantum devices like quantum gates, filters and 
qubits when the boundary parameters become tunable and controllable [T5] . 



9 Conclusion 

We have seen how the ambiguity in the self-adjointness of a Schrodinger Hamiltonian with 
a singular potential can be understood, and how it can be technically treated in a framework 
that requires only a reduced amount of calculational efforts. 

Naturally, what has been presented here is only a narrow and practice-oriented extract from 
the extensive research field of boundary conditions and singular potentials and, more generally, 
of self-adjoint extensions. For further reading, in addition to the works already referred to, one 
can start with consulting, for example, [171 \18\ \19\ 1201 \21\ [22] and references therein. 

To close the discussion with the message: the self-adjointness ambiguity caused by irregulari- 
ties of the potential can be interpreted such that the system carries extra physical properties 
which can not be expressed via the potential function but through a boundary condition at each 
irregularity. These properties are to be fixed either by some additional theoretical knowledge 
or physical information about the system, or by measurement, fitting the unknown boundary 
parameters to experimental data. 

12 Writing -% as Kl+B [ ambiguity in the self-adjointness occurs for I < \- In parallel, the Hamiltonian will 
be bounded from below when g > — ^ (when I is real). 
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